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Abstract 

By means of adding a collision process between the ball and racket in double pendulum model, 
we analyzed the tennis stroke. It is possible that the speed of the rebound ball does not simply 
depend on the angular velocity of the racket, and higher angular velocity sometimes gives lower 
ball speed. We numerically showed that the proper time lagged racket rotation increases the speed 
of the rebound ball by 20%. We also showed that the elbow should move in order to add the 
angular velocity of the racket. 
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I. INTRODUCTION 


The double pendulum has been studied as an example of the chaotic motion in physics [I]- 
[Hj. If we consider the first cycle of the pendulum, the double pendulum model holds appli¬ 
cation in sports such as golfj2-[7], baseball [8], and tennis [9]. Since the double pendulum is 
not a simple linear system, the motion of the pendulum can not be optimized as a simple 
analytic form. The swing pattern utilised to maximize the angular velocity of the hitting 
rod such as racket, bat, and club has been analyzed on the assumption that the angular 
velocity is the dominant factor for the speed of the rebound ball[S]. 

For a tennis stroke, the two rods for double pendulum are an arm and a forearm for first 
rod and a racket for the second rod. In our model, we added the collision process between 
the ball and the racket. Without the collision process j9], there is no criteria to attain high 
speed of the rebound ball except for the angular speed of the second rod, racket. If we set 
the impact angle of the first rod at which the ball hits the racket, the speed of the rebound 
ball is mainly dependent on the angular velocity of the second rod. On the other hand, if we 
release the impact angle of the first rod, the speed of the rebound ball does not remain as 
a simple function of the angular speed of the hitting rod. We showed that the speed of the 
rebound ball is different even though the angular velocities of the second rod at the contact 
time are same. Furthermore, considering the whole stroke, the maximum angular velocity 
for the lower speed of the rebound ball is greater than the maximum angular velocity for 
the higher speed of the rebound ball. Therefore, to get maximum speed of the rebound ball, 
it’s not sufficient to set the condition to generate high angular velocity of the hitting racket. 

The collision between the racket and ball has been studied in various ways In our 

simple collision model, we assumed that the racket is a simple one dimensional rod without 
any nodal motion and the collision occurs in one dimension. This assumption is valid if the 
racket and a ball moves in the same line for the short collision time. Although, our model 
does not give any detailed information on the collision about the effect of the tension of the 
string and the mass distribution of the racket, however provides some insights on the proper 
swinging pattern to get maximum speed of the rebound ball 

In this article, we also analyzed the tame lagged torque effect for the double pendulum 
system. By applying time independent constant torques on first rod and the second rod, 
the speed of the rebound ball can be calculated for certain initial conditions. In the same 
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condition, if we simply hold the racket for a short time without enforcing a torque and with 
subsequent application of torque, the naive intuition estimates decrease in the speed of the 
rebound ball. The reason behind this is application of less energy for the double pendulum 
system. However, the speed of the rebound ball increase by 20% by choosing the proper 
delay time. It’s mainly because the double pendulum system is not a simple linear system. 
Adding energy to the double pendulum system does not directly increase the speed of the 
rebound ball. We also analyzed movement of the elbow to which the first rod of the double 
pendulum is attached. At a first glance, if we add extra movement towards the rebound 
ball’s direction at contact time, the speed of the ball increases; however, with varying results. 
Apparently, it becomes clear that, the double pendulum system is really a nonlinear system. 

The present paper is organized as follows: In Section II, we introduce a double pendulum 
system including the collision process. In Section III, the differential equations obtained 
in section II are solved, we numerically showed that the speed of the rebound ball does 
not simply depend on the angular speed of the racket. For some cases, the higher angular 
velocity gives lower speed of the rebound ball. In Section IV, we analyzed the dependence of 
the racket mass and length of the first rod. The general properties of swing system has been 
demonstrated. When we applied a time dependent torque on the racket, we could increase 
the speed of the rebound ball. The time lagged rotation of the racket was analyzed in section 
V. In section VI, the elbow movement is analyzed to add additional speed to the ball. In 
Section VII, we summarize the main results and discuss the application of our results. 

II. RACKET SYSTEM INCLUDING THE COLLISION WITH A BALL 

The geometry of the double pendulum model for the swing of a racket is shown in Fig 
0( a). Though, this geometry is for the left-handed player if we see from the +z direction, 
it’s originally related to the real double pendulum problem in the gravitational field [8j. 
Our basic model and some notations are closely related with those in work reported by Rod 
Cross in]. The elbow moves in the xy plane and the arm and the racket also moves in xy 
plane. We also modeled the arm and the forearm as a simple uniform rod with mass Mi and 
length Li. The racket including the hands is also treated as a uniform rod with mass M 2 
and length L 2 . The arm and the racket rotate in a clockwise direction in a plane at angular 
velocities oji = —dO/dt and U 2 = —d(f)/dt, respectively. If we assume that the velocity of the 
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elbow is (14®, V^) then the velocities of the center of the first rod (arm and forearm system) 
and the second rod (hand and racket system) becomes 

Vxi = Vx ~ h\Ui cos 6 
Vyi = Vy — hiUJi sin 9 
Vx2 = Vf — LiUi cos 9 — h2U!2 cos 0 

Vy2 = — LiUi sin 9 — h2U)2 sin 0 ( 1 ) 

,where {vxi,Vyi) , {vx2,Vy2) are the velocity of the center of the mass of the first rod and 

the second rod, respectively. The center of masses are located in the middle of the rod, 

hi = Li /2 and ^2 = T2/2, since we assumed uniform rods. 



FIG. 1 : Schematic diagram of arm and racket system, (a) Length of the first rod (forearm and 
arm system) is Li, and the length of the racket is L2. The angles 0 and </> are defined from the 
—y direction, (b) The Force on two rods. {Fxi,Fyi) is the force acted on the first rod at the joint 
between elbow and the first rod. {Fx2, Fy2) is the force acted on the second rod at the joint between 
the two rods. F^oi is the force on the second rod by a ball. 

Let the force from the elbow to the first rod on joint point between the elbow and the 
first rod be {Fxi,Fyi)^ and the force from the first rod on joint point between two rods be 
{Fx2i Fy2), then the equations of the motion for two center of masses become 

Mi^ = Fxi - Fx2 
at 

( 2 ) 
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,, dVx 2 ^ ^ 

M2 — 7 — = Fx 2 + Fcol 


dt 


Mo 


dv. 


J/2 


dt 


= F 


y2 


( 3 ) 


The two forces {—Fx2, —Fy2) are reaction force of {Fx2, Fy2) from the second rod. We added 
the force Fh from the ball to the second rod at the point of hn from the center of the mass. 
This force from the collision between the racket and the ball changes the torque equation in 
[ 9 ] and the two torque equations becomes 
di 02 


L. 


2,cm 


^ l,cm' 


dt 

duji 

dt 


= C2 + Fx2h2 cos (p + Fy2h2 sin p - Fcoihn cos p 
= Cl — C2 + Fxihi cos 6 

T Fyihi sin 0 T Fx2(^Fi — hi^ cos 0 T Fy2(^Li — sin 0 


( 4 ) 


( 5 ) 


,where Ci is the torque on joint point between the elbow and the hrst rod and C2 is the 
torque on the joint point between the two rods. From the Eqs. i-i we obtain two equations 
for the time derivative of two angular velocities as follows 

duJi P{l2,cm + M2/I2) - (3M2/I2T1 cos(0 - 0) + S'! 


dt ~ {l2,cm + M2hi)(/i,,^ + Mihl + M 2 LI) - MihlLl cos^iP - 6) 
dc02 _ Q{Ii,cm + Mihl + M2Ll)-PM2h2Licos{P-9) + S2 
dt ~ {l2,cm + M2hl){Ii^,m + Mihj + M 2 LI) - Mlh^L? cos^{p - 6) 


( 6 ) 


where tty = and 

P — Cl — C 2 — M2h2LiUjlsm{p — 6) + (Mihi + M2Li)[axCos9 + ttySmO] 

Q = C 2 + M2h2LiU)lsm{p — 9) + M2h2[axCosp + aysinp] (7) 

P and Q are similar to the results in [9j with setting g = 0, and the collision force from the 
ball add following two terms 

51 = FcoiLi[{2l2,cm +h 2 {h 2 -hn)M 2 ) cos 9 - h 2 ih 2 +hn)M 2 cos{9 - 2p)] 

5 2 = Fcoi[{2hH{Ii,cm F Mihl + M2L‘1) 

+ /i 2 ( 2 /i^cm + 2Mihl + M 2 L\)) cos P - h 2 M 2 L\ cos(20 - p)] (8) 

The collision force defines the motion of the ball as follow , 
d‘^Xb „ 

7.0 P^col 


dP 


— —fki^H — Xb) 


xh = xe F Li sin 9 + (h2 F hn) sin p 
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( 9 ) 

( 10 ) 











,where xh is the x component of the hitting part of the racket which locates h2 + hn from 
the bottom of the racket, and xe is the x component of the joint between the elbow and 
arm. 

The force between the racket and the ball should be repulsive, so that the force form 
becomes 

fkixH - Xb) = k{xH - Xb) if Xh > xt 

= 0 if Xb> Xh- ( 11 ) 

If we consider this harmonic force between the ball and the second rod, the period of the 
oscillation is 

T= =. (12) 

/ k(mb+M2) 

Y rrn,M2 

Since the half of this period is the collision duration between the ball and the racket, we 
controlled the k value from 188 A^/m to 19000 A^/m. Subsequently, the collision time varies 
from 50 ms to 5 ms. However, if the contraction length is large, the Hook’s model is not valid 
for the ball and the racket system. Then the force can be rewritten as kxmsin(^ ) ra . or 
tan(^^) in order to limit the maximum contraction length. However, in our numerical 
calculation we restricted our system in order to follow the Hook’s rule. 

In this article, we assumed a simple model for the ball and racket collision. In our model, 
the ball is assumed to hit the racket when the racket is parallel to y axis (0 = 0 ). In addition 
to this, the ball is assumed to be moving in x axis. In this case, we numerically calculated 
0 {t) and fit) till the time when the racket is parallel to the y axis ( 0 (to) = 0 ). With 
these numerical results, we set new initial conditions just before the collision. We assumed 
that the collision occurs aX t — Iq — T/ 2 ^ where T is the period of the harmonic oscillator 
system between the racket and the ball. At this time, we numerically solved the ball and 
two rod system attached to the moving elbow numerically till the ball is rebound from the 
racket. We determined the speed of the ball at that moment. Based on our simple model, 
the coupling constant k determines the collision duration, but the speed of the rebound ball 
is not altered very much. The reason for such a behavior is consideration of one dimensional 
collision. However, the purpose of this article was to find the optimum path for the swing of 
the racket, we did not extend our model to two dimension. We only restricted our numerical 
conditions to render presence of our setups in the valid region. 
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III. DIFFERENT BALL SPEED WITH THE SAME ANGULAR VELOCITY OF 


THE RACKET 

We assumed that an arm and a forearm forms a simple rod and the moment of inertia 
about the center of mass is /i,cm = MiL\/2. The moment of inertia about the center of mass 
of the hand-racket system is also assumed as l 2 ,cm = M 2 L 2 / 2 . Although this assumptions 
are not enough to study the tennis stroke in detail, we focused our attention on the double 
pendulum model, which has its application in the tennis stroke. In this section, the torques 
applied to the first forearm system Ci and to the racket system C 2 were set by 25V and 
2 .5V as described in Rod Cross work [9]. The velocity {m/s) of the elbow is assumed as 
follows [9] , 

(t) = —33t + 69t^ 

v^{t) = -A2t + 174^2. (13) 

In Fig. we plotted the speed of the rebound ball from the racket as a fnnction of the 
initial angles (6*0, (5o). The angle S = (p — 0 is defined as the angle between the first rod and 
the racket, ho is the initial angle at t = 0. We assumed that the length of the first rod 
(an arm and a forearm system) as Li = 0.3m and the mass of the first rod as Mi = 2.0kg. 
And the length of the racket system as L 2 = 0.7 and the mass of the racket system as 
M 2 = 0.3kg. We assnmed that the initial velocity of the ball just before the contact is 5m/s 
to the positive x direction. The xh , a hitting position from the center of mass assumed 
to be as 0.15m. The length of the first rod (an arm and forearm system) assumed to be as 
0.3m. This length is projection length in x — y planes. We plotted initial conditions which 
gives maximum speed of the rebound ball. As 9o changes from 60° to 120°, the becomes 
smaller and reaches 20°. In Fig. we also plotted the speed of the rebound ball from the 
racket as a fnnction of the initial angles (6*0, (jo)- We only changed the length of the first rod, 
in other words we extended the distance between the elbow and the hand. At this time the 
(jo is smaller as compared to the case when Li = 0.3m. When the initial angle Oq is 120°, the 
racket and the first rod should be in the same line so as to get the maximum speed of the 
ball. In our model, the speed of the rebound ball is not a simple function of the angular 
velocity 002 . In Fig. we plotted the angular velocities ui and UJ 2 for two cases. For the case 
a) the initial angle (6*0, (jo) is (120.0°, 81.5°) and the contact time ta = 0.198s and the angnlar 
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FIG. 2: The speed of the ball in km/h unit as a function of the initial angle of the racket and arm 
with respect to the —y axis. Oq : initial angle of the arm. 5o ^ initial angle between the racket and 
the first rod. The length of the first rod is 0.3m. Red line indicates two initial angles which gives 
the maximum speed. 



FIG. 3: The speed of the ball in km/h unit as a function of the initial angle of the racket and arm 
with respect to the —y axis. Oq : initial angle of the arm. Sq : initial angle between the racket and 
the first rod. The length of the first rod is 0.4m. Red line indicates two initial angles which gives 
the maximum speed 
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velocity is coi = 11.34rad/s. For the case b), the initial angle (0o,(^o) is (75.0°,39.0°) and 
the contact time tf, = 0.220s and the angular velocity uji = 17.9rad/s. However, the angular 
velocities u )2 at the contact times have the similar value as 37.98rad/s. Furthermore, the 
angular velocity U 2 b increases even after the collision time. The interesting thing is that the 
speed of the rebound ball for the case b) is smaller than the speed of the ball for the case 
a). The speeds are IbQkm/h and 141/cm/h, respectively. 

If we only check the angular velocity a; 2 , we may conclude that the case b) gives higher 
speed of the rebound ball. But the angular velocity is not the entire factor to determine the 
speed of the ball. This can be explained if we examine the angle (dc) of the first rod when 
the racket contacts the ball. are 11.5° and 87.1° for the cases a) and b), respectively. 
The speed of the ball {vout) is also a functions of 9c as well as ui and u) 2 . In order to get the 
maximum speed of the rebound ball, the double pendulum system should be examined as a 
whole system. 



FIG. 4: The angular velocities loi and L 02 for two cases. The initial angle (0O)<io) is (120.0°, 81.5°) 
for the case a) and (75.0°, 39.0°) for the case b). ta and 4 are the contact time. 


IV. DEPENDENCE OF THE RACKET MASS AND THE DISTANCE BETWEEN 
THE ELBOW AND THE HAND 

The mass of the racket-hand system is assumed to be 0.3kg, but the mass of the racket 
may be changed. All the movements are assumed to be in the xy plane in our model, the 
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projected length between the elbow and the hand in xy plane can be changed by controlling 
the angle between the forearm and arm in actual tennis stroke. In Fig. we plotted the 
speed of the rebound ball from the racket as a function of the racket mass (M 2 ) and the 
projected length in xy plane between the elbow and the hand Li. If the mass of the racket is 
about 0.2%, the speed of the rebound ball is decreased as the projected length between the 
elbow and the hand is increased. However, if the racket mass is getting heavier, the speed 
of the rebound ball increases. Since an actual mass of the tennis racket is around 0.3%, and 
L\ is limited, the speed of the ball is restricted. 

In Fig. 1^ we plotted the angle {6c) of the forearm system at the impact time as a function 
of the racket mass (M 2 ) and the projected length in xy plane between the elbow and the 
hand Li Regardless of the racket mass (M 2 ), the forearm angle 6c decrease to 20° as the 
length Li increase. This results shows that for the player with folded arm whose effective 
projected length of the forearm system is small, the contact angle (do) should be about 90°. 
If the length L\ and the contact angle do are reduced the ball should be hit relatively close 
to the body. This tendency may explain the impact point of the ball in tennis stroke and 
golf. 



FIG. 5: The speed of the rebound ball from the racket as a function of the racket mass (M 2 ) and 
the projected length in xy plane between the elbow and the hand Li. The unit of Vgut is km/h 
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FIG. 6: The angle of the forearm system at the impact time as a function of the racket mass (M 2 ) 
and the projected length m. x — y plane between the elbow and the hand L\. 

V. TIME LAGGED TORQUE VS. CONSTANT TORQUE. 

In actual tennis stroke, most of the players use time lagged racket movement. They 
intentionally keep the racket back as forearm rotates then they start to move the racket to 
get a high angular velocity UJ 2 . The main difference of double pendulum when compared to 
the single pendulum, is the separate movement of the first rod (arm and forearm system) 
and the second rod (the racket system). In Fig. we plotted the speed of the rebound ball 
as a function of increase in the torque C 2 and the delay time r. We set the projected length 
in xy plane between the elbow and the hand Li as 0.4m. The torque applied to the racket 
at the joint of two rods, varies from 0 to ION. The time delay r is the starting time at 
which we applied the torque to the racket. After waiting for r second, the torque suddenly 
changes from zero to a certain value C 2 till 0.5ms in our numerical calculation. This time is 
far behind the contact time (around 0.2ms). At first glance, if we start to apply the torque 
early, the angular velocity of the racket may accelerate little bit more. It’s simply because 
the acceleration of the angular velocity depends on the torque. If the torque C 2 is less than 
1.56A, the maximum speed of the rebound ball is obtained when the delay time r is zero 
as we expected. However, for the higher torque C 2 > 1.56A,the numerical results in Fig. 
are different from our simple intuition. For a given torque value C 2 > 1.56A, there always 
exist a certain time delay r at which the speed of the rebound ball is maximum and the r 
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is not zero. In other words, the important thing to get high speed of the rebound is not the 
total amount of impulse (torque x (applied time)), but the timing when the torque starts. 


'^out 



FIG. 7: The speed of the rebound ball as a function of the torque on the racket {C 2 ) and the 
time lag r. The initial angle (5o is 26.1°, the which angle gives the maximum speed for the angle 
(00 = 90°). Li = 0.4m 


If the torque C 2 is lagged, the initial angle 5 = 26.1° to get a maximum speed of the 
rebound ball for the fixed initial angle Oq = 90° is no more the optimum condition. In order 
to see the time lag effect clearly, we set the torque C 2 = . In Fig. we plotted the speed 

of the rebound ball as a function of the time lag r and initial angle ho with an initial angle 
9q = 90°. The new optimum condition to get maximum speed of the rebound ball is the 
time lag r = 0.185s and the initial angle ho = 52.5°. In Fig. we plotted the speed of the 
rebound ball as a simple function of time delay r. If the time lag r of the torque 6*2 is r/, = 0, 
the rebound ball speed is vl = 136.6A:m/h. And Vout has its maximum [vh = IQAAkm/h) at 
th = 0.185s. The speed increases by about 20% by adjusting the time delay with the same 
magnitude of the torque € 2 - The angular velocity of the forearm system and the racket 


system is shown in Fig. 10, For two time delay tl and th, the difference in the contact time 


is less then 10ms. The dashed line in Fig. [ 10 ] indicates the time at which the ball collides 
with the racket. We note that the racket hit the ball before the racket has its maximum 
angular velocity. In earlier work jS] , the authors have analyzed the double pendulum system 
in order to get the maximum angular velocity u) 2 . However, this may mislead the double 
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FIG. 8 : The speed of the rebound ball (km/h) as a function of the time lag r and initial angle (5o 
with an initial angle Oq = 90° 

pendulum system. Considering Fig. it is clear that higher angular velocity gives lower 
speed of the rebound ball. 


In Fig. when the angular velocity U 2 h has its maximum, the angular velocity tui is 
almost zero. This is explained for the double pendulum model of tennis stroke. The angular 
momentum and energy of the first rod (arm and forearm system) were totally transferred 
to the racket system in order to get maximum angular velocity of the racket. 

We plotted the first rod (forearm and arm system) and the racket system for the time 
delays represented as th and tl in Fig. and [T^ respectively. When the first rod starts 


to rotate, the racket stayed back in Fig. till time th- After applying the torque 6 * 2 , 
the racket suddenly starts to rotate, and hits the ball with the angle 6c = 53.3°. The 
racket and the first rod are shown in red at the contact time. Comparing the stroke with a 
constant torque (Fig. [^, the racket rotates more rapidly at the contact time. In double 
pendulum model for the stroke, this phenomena was expected qualitatively. In our model 
we quantitatively demonstrated that why time lagged stroke is needed and the extent to 
which the velocity can be increased. , 
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FIG. 9: The speed of the rebound ball as a simple function of time delay r. At r = 0.185s, the 
Vout has it’s maximum . 



FIG. 10: The angular velocity of the forearm system uii and the racket system {omega 2 ). 

{uiil^^2l) indicates the angular velocities when the time delay is tl, and {uiihi^2h) indicates 
the angular velocities when the time delay is th 

VI. NEED MORE IN ORDER TO INCREASE THE SPEED OF THE REBOUND 
BALL. 

In the previous section, we noted that the important thing is not the total impulse to 
accelerate the racket, but the timing in order to increase the speed of the rebound ball. 
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FIG. 11: Time dependent trace of the tennis racket system (black bar) and the arm and forearm 
system (blue bar), when the time delay is th- The motions are captured from t = 0 to t = 0.325 
evenly. The red one is the racket position at the contact time. 



FIG. 12: Time dependent trace of the tennis racket system (black bar) and the arm and forearm 
system (blue bar), when the time delay is tl = 0. The motions are captured from t = 0 to t = 0.325 
evenly. The red one is the racket position at the contact time. 

Now, we analyze the influence of the elbow movement on the speed of the rebound ball. We 
made a simple model to add an additional movement of the elbow. The added force around 
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the time th as follows, 

F^{t) = Fo cos -^F exp(-(t - {th + rF)f/tl) 

Fy{t) = Fo sin exp(-(t - {th + tf)Y /t\) (14) 

where 4/^? is the angle measured from the +x axis counterclockwise, we also set Fq = 20N, 
and the width = 0.05s. Since the ball is rebound to —x axis, the proper direction of 
the force seems to be in the —x direction. We plotted the speed of the rebound ball when 
the extra force is applied in the direction of d'i;’ and with extra time delay rp in Fig. 
Numerical results show that we can obtain the maximum speed of the rebound ball when 
the angle is 76°. When the direction of the force 'i/p = 180°, the speed of the rebound ball 
is almost same to the speed of the ball without the additional force. In other words, the 
external movement towards —x direction does not give any additional speed to the rebound 
ball. In Fig. we plotted the trajectory of the elbow from the time t = 0 till the contact 
time. F/vf is the path when no additional force is added. Pmax, Pmin are elbow’s trajectory 

when the direction of the force ^/f are 76° and, 275°, respectively. Table |I] demonstrates 

the Vout and the angular velocity U 2 for 4-cases. The movement with 4/^ = 76° is slightly 
backward and mainly perpendicular to the direction of the rebound ball. From this result, 
we conclude that the main factor for the increased speed of the rebound ball is not the linear 
momentum added to the racket but the angular momentum added to the racket system. 



FIG. 13: The speed of the rebound ball as a functions of the extra force with a direction of ^^f 
and with a extra time delay Tp 
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FIG. 14: The trajectory of the elbow for four cases from the time t = 0 till the contact time. Pnf 
is the path when no additional force is added. Pforward-, Pmax-, Pmin are elbow’s trajectory when 
the direction of the force are 180°, 76°, 275°, respectively. 


TABLE I: Additional Force dependent output 


Fo (N) 

^^(°) 

Vout (km/h) 

(jj 2 rad/s 

20 

180 

163.3 

39.1 

20 

275 

158.8 

39.2 

20 

76 

170.6 

52.2 

0 

0 

164.4 

44.7 


VII. CONCLUSION AND DISCUSSION. 

The double pendulum model is applied to the baseball, tennis, and golf. It analyzes the 
swing pattern to maximize the angular velocity of the hitting rod such as racket, bat, and 
club on the assumption that the angular velocity is the dominant factor to attain speed of 
the rebound ball. If we set, Oc, the angle of the first rod at the impact time, it seems to 
be obviously reasonable. On the other hands, if we release the angle 6c, the speed of the 
rebound ball is not a simple function of the angular speed of the hitting rod. The speed of 
the rebound ball is different even though the angular velocities at the contact time are same 
in Fig. 1^. Furthermore, considering the whole stroke, the maximum angular velocity 002 for 
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the lower speed of the rebound ball is greater than the maximum angular velocity (j 02 for the 
higher speed of the rebound ball. Therefore, to attain maximum speed of the rebound ball, 
it’s not sufficient to set the condition to generate high angular velocity uj 2 of the hitting rod. 

In the double pendulum, the efficient way to generate high angular velocity of the hitting 
rod is the energy transfer from the first rod to the second hitting rod. In other words, when 
the angular velocity (Ui of the forearm system is zero, the angular velocity of the hitting 
racket has maximum value as shown in Fig. [l^ We analyzed the time lagged torque effect 
for the double pendulum system. With applying constant torques C\^C 2 on the forearm 
system and on the racket, respectively, the speed of the rebound ball can be calculated. For 
the same condition, if we simply hold the racket for a short time without enforcing a torque 
C 2 , then applying the torque C 2 at the proper time {t = th), the speed of the rebound 
ball increases by 20% as can be seen in Fig. The reason is mainly because the double 
pendulum system is not a simple linear system. Adding the energy to the double pendulum 
system does not directly increase the speed of the rebound ball. 

We also analyzed the elbow movement effect. In addition to the velocity of the elbow for 
the medium pace forehand, we added extra movement of the elbow. At a first glance, if we 
add extra movement towards the rebound ball’s direction, the speed of the ball is increased. 
But the double pendulum system is not a simple linear system. When the direction of the 
elbow movement is perpendicular to the ball’s direction, the speed of the rebound increases. 
Actually the direction is towards the center of the elbow’s circular movement. In other 
words, the added centripetal force does not add not the linear momentum of the racket, but 
the angular velocity of the racket. 

Although our collision model is applied in one dimension, this collision process allows us 
to analyze the double pendulum system in a more realistic manner. We showed that the 
speed of the rebound ball does not simply depend on the angular velocity of the racket. 
The increase in the ball speed by the proper time lagged racket rotation was numerically 
studied. The elbow movement for adding the ball’s speed was counter intuitive. The addition 
of simple linear momentum to the elbow is not important; however, the elbow should move 
in order to add angular velocity to the racket. 

In actual tennis stroke, the motion occurs in three dimensions and the magnitudes of the 
forces and torques are dependent on the muscle shape and movement. We did not include 
any bio-mechanical information such as pronation and we did not include the spin of the 
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ball in any way. The numerical data may also be not suitable for some players. However, 
our study on the double pendulum system for tennis stroke provides some insights to attain 
an efficient way to stroke a tennis ball. 
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